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ABSTRACT 

We report on optical integral field spectroscopy of two unrelated blue compact 
galaxies mapped with the 13 x 13 arcsec 2 VIMOS integral field unit at a resolution 
of 0.33 x 0.33 arcsec 2 . Continuum and background subtracted emission line maps in 
the light of [O in] A5007, Ha, and [N n] A6584 are presented. Both galaxies display 
signs of ongoing perturbation and/or interaction. UM 420 is resolved for the first 
time to be a merging system composed of two starbursting components with an 'arm- 
like' structure associated with the largest component. UM 462 which is a disrupted 
system of irregular morphology is resolved into at least four starbursting regions. 
Maps of the Ha radial velocity and FWHM are discussed. No underlying broad line 
region was detected from either galaxy as the emission lines are well-fitted with single 
Gaussian profiles only. Electron temperatures and densities as well as the abundances 
of helium, oxygen, nitrogen, and sulphur were computed from spectra integrated over 
the whole galaxies and for each area of recent star formation. Maps of the O/H ratio 
are presented: these galaxies show oxygen abundances that are ~ 20 per cent solar. No 
evidence of substantial abundance variations across the galaxies that would point to 
significant nitrogen or oxygen self-enrichment is found (< 0.2 dex limit). Contrary to 
previous observations, this analysis does not support the classification of these BCGs 
as Wolf-Rayet galaxies as the characteristic broad emission line features have not been 
detected in our spectra. Baldwin-Phillips- Terlevich emission line ratio diagrams which 
were constructed on a pixel by pixel basis indicate that the optical spectra of these 
systems are predominantly excited by stellar photoionization. 

Key words: galaxies: abundances - galaxies: individual (UM 420, UM 462)- galaxies: 
kinematics and dynamics - galaxies: starburst - galaxies: dwarf 



1 INTRODUCTION 

Blue compact dwarf galaxies (BCGs) offer a means of ex- 
ploring star^onnatic4i_Jji low-m ass, low-metallicity (1/50- 



1/3 Zq Kunth & Ostlin (20001) systems and, by analogy, 
the chemically un-evolved systems in the high-z primordial 
universe. Typically, they are characterised by their strong 
emission lines superimposed on a faint, blue continuum, re- 
sulting from one or more recent bursts of star-formation. 
These compact, gas-rich objects are estimated to ha ve high 
star- formation rates in the order of 0.1 to 1 M fl yr _1 (Fanelli 



et al.|l988| with typical H I masses of Ma; ~ 10 7 - 10° M 
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( |Thuan fc Martin|1981[ ), and are ideal laboratories for many 
topics related to star-formation, partly because they do 
not display more complicated phenomena, such as density 
waves, that operate in larger galaxies ( [Cairos et al.|[200T| ). 
This is a study of two BCGs, UM 420 and UM 462, first 
discovered by the University of Michigan survey for extra- 
galactic emission-line objects ( |MacAlpine fc Lewis|[l978"| ). 
Their individual properties may provide insights into the 
role of galaxy interactions in chemical recycling and star- 
bursting episodes. Several studies have suggested that inter- 
actions with other systems are a contributing factor to large- 



scale starbursts in dwarf galaxies l 


Brinks 


1990| Mendez 


& Estcban 2000. Iglcsias-Paramo & Vflchez 


|2001| |Verdes- 


Montenegro et al.| 


2002| |Tran et al.||2003 


1. Hence it may be 



no coincidence that both galaxies studied here display signs 
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of interaction and/or perturbation in the form of tails, mul- 
tiple nuclei, or disrupted morphology. 

UM 420 was reported by |L6p cz-Sanchcz & Estcban 
(2008) as having an irregular and elongated morphology, 
with two long regions extending in different directions 
from the brightest central region. A bright spiral galaxy, 
UGC 1809, lies 16.5" west of UM 420, but at less than half 
the redshift of UM 420 is not physically associated to it and 
cannot cause an effect. UM 420 was once thought to be a 
large H II region near the edge of UGC 1809's spiral arms, 
and as a result the SIMBAD data base currently incorrectly 
lists its RA and DEC as those of UCG 1809 (the correct RA 
and D EC for UM 420 are given in Table [T| . |Pustilnik et aT 



( 2004 1 describe UM 420 as being a possible 'sinking merger', 



i.e. harbouring a satellite galaxy that is sinking into a larger, 
gas-rich companion. 

UM 462 with UM 461 form part of an interact i ng bi- 
nary pair separate d by only ~50" ( Taylor et al.|1995 |Telles 



& Terlevich 19951. In a Hi survey of H II galaxies J with 
companions ( [Taylor et al~1p95| , UM 461 was observed to 
have a tidal arm extending in the direction of UM 462 and 
outer Hi contours that are clearly disturbed by its nearby 
neighbour. |Telles fc Terlevich| ( |1995[ ) interpreted UM 462 as 
being part of a group of H II galaxies, UM 461, UM 463 and 
UM 465, whose star-bursting episodes are synchronised on 
time scales of less than 10 7 years, despite being separated by 
~l-2 Mpc. Following the classification system of |Loose fc| 
Thuan|([T985l),|Cair6s et al.|(|2001[) classified UM 462 as being 



of 'iE' type, i.e. showing a complex inner structure with sev- 
eral star-forming regions superimposed on an extended regu- 
lar envelope. This type of structure is commonly interpreted 
as a sign of interaction. Spectroscopically, both galaxies have 
been classified as Wolf-Rayet (WR) galaxies (this is however 
disputed in the present work; see Sec. 5.3), with broad He II 
A4686 being detected in their spect ra (|Izotov fc Thuan|1998| 
Guseva et al.|2000 l. Guseva et al. (20001 also claimed to see 
broad C IV A4658 and C IV A5808 emission in the spectra 
of UM 420, indicating the presence of WCE stars, although 
the detections were doubted by Schaerer et al. ( 1999 1 on 



the basis of the S/N ratio of the Guseva et al. (20001 spec- 
trum. The relative numbers of WR stars to O-type stars 



were estimated by Guseva et al. (20001 as 0.03 and 0.005 



for UM 420 and UM 462, respectively. Schaerer et al. (1999) 



have defined WR galaxies as being those whose ongoing or 
recent star formation has produced stars sufficiently mas- 
sive to evolve to the WR stage. However, the same authors 
also note that the definition of a WR galaxy is dependent 
on the quality of the spectrum, and the location and size of 
the aperture. 

UM 420 falls within a small sub-set of BCGs reported by 



Pustilnik et al. (20041 as having anomalously high N/O ra- 



tios (^ 

licity (12+log(Q/H)=7.93; |Izotov , 

after). They attributed this over- abundance of nitrogen to 
the effects of N-rich winds from large numbers of Wolf-Rayet 



0.5 dex) compared to BCGs of similar oxygen metal- 

; Thuan|1998] IT98 here- 



(WNL-type) stars. A previous IFU study by James et al 



(2009) (J09 hereafter) of another BCG in this 'high N/O' 



H II galaxies are galaxies with spectra simila r to H II regions 



but more metal-rich than blue compact dwarfs (Kunth & Ostlin 



2000) 



group, Mrk 996, derived a significantly higher O/H ratio 
than previously and an N/O ratio for its narrow-line com- 
ponent which is now typical for its metallicity. Mrk 996, how- 
ever, was shown by J09 to also contain an extremely dense, 
extended broad-line component displaying a N/O ratio up 
to 20 times higher than the galaxy's narrow-line component; 
this large N/H and N/O enrichment was attributed to the 
cumulative effects of a large population of ~2600 WNL-type 
stars whose number was estimated from fitting the broad 
emission features at 4640 and 4650-A with theoretical WR 
spectral templates. By conducting a similar study of UM 420 
we aimed to investigate its N/O ratio status, and to study 
potential effects to the gas abundance patterns due to WR 
stars. On the other hand, UM 462 provides a useful compar- 
ison, as it has been reported to have a rather normal N/O 
ratio close to the average observed for many BCGs at its 
published metallicity of 12 + log(0/H) = 7.95 (IT98). 

In order to understand interaction-induced starbursts, 
it is essential to gain spatial information regarding the prop- 
erties (e.g. physical conditions, chemical abundances) of the 
starburst regions and any spatial correlations that may hold. 
This type of information has been limited in previous long- 
slit spectroscopic studies of H II galaxies. In this paper we 
present high resolution optical observations obtained with 
the VIMOS integral field unit (IFU) spectrograph on the 
ESO 8.2m Very Large Telescope UT3/Melipal. The data 
afford us new spatiokinematic '3-D' views of UM 420 and 
UM 462. The spatial and spectral resolutions achieved (Sec- 
tion 2.1 1 allow us to undertake a chemical and kinematical 



analysis of both systems, providing a clearer picture of the 
ionised gas within their star-forming regions. We adopt dis- 
tances of 238 Mpc and 14.4 Mpc for UM 420 and UM 462, re- 
spectively, corresponding to their redshifted velocities mea- 
sured here of +17,454 and +1057 km s _1 (z = 0.060604 and 
0.003527, as given in Table fl|, for a Hubble constant of H 
= 73.5 km s" 1 Mpc" 1 ( jDeBernardis et al.|2008 |. 



2 VIMOS IFU OBSERVATIONS AND DATA 
REDUCTION 

2.1 Observations 

Two data sets for each galaxy were obtained with the Visible 
Multi-Object Spectrograph (VIMOS) IFU mounted on the 
8.2m VLT at the Paranal Observatory in Chile. All data sets 
were taken with the high-resolution and high-magnification 
settings of the IFU, which resulted in a field-of-view (FoV) 
of 13" x 13" covered by 1600 spaxels (spatial pixels) at a 
sampling of 0.33" per spaxel. The data consist of high- 
resolution blue grism spectra (HRblue, ~0.51 A/pixel, 2.01+ 
0.23 A FWHM resolution) covering 4150-6200 A and high- 
resolution orange grism spectra (HRorange, ~0. 6 A/pixel, 
1.92 ± 0.09 A FWHM resolution) covering 5250-7400 A. 
Spectra from both grisms for each galaxy are shown in Fig- 
ure [l] The observing log can be found in Table [2] Four ex- 
posures were taken per grism, per observation ID. The third 
exposure within each set was dithered by +0.25" in RA and 
by +0.61" in DEC for both galaxies (corresponding to a 2 
spaxel offset in the X-direction) in order to remove any bro- 
ken fibres when averaging exposures. All observations were 
taken at a position angle of 20° for UM 420 and 22° for 
UM 462. 
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Table 1. General Parameters of UM 420 and UM 462 





Coordinates (J2000) 








Name 


a 8 


z a 


Distance (Mpc) 


Other Names 


UM 420 


02 21 36.2 00 29 34.3 


0.060604 


238 




UM 462 


11 53 18.5 -02 32 20.6 


0.003527 


14.4 


Mrk 1307, UGC 06850 



a Derived from the present observations 
Table 2. VIMOS IFU observing log 



Observation ID 


Date 


Grism 


Exp. time (s) 


Airmass range 


FWHM seeing (arcsec) 








UM 420 






250536-9 


14/11/2006 


HR blue 


4 x 402 


1.11-1.40 


0.85 


250534,250535 


14/11/2006 


HR orange 


4 x 372 


1.13-1.36 


0.79 








UM 462 






250563-6 


22-23/03/2007 


HR blue 


4 X 402 


1.60 - 1.10 


1.02 


250559,250562 


22-23/03/2007 


HR orange 


4 x 372 


2.16 - 1.73 


0.74 



2.2 Data Reduction 

Data reduction was carried out using the ESO pipeline 
via the GASGANcj^] software package and followed the se- 
quence outlined by J09. The reduction involves three main 
tasks: vmbias, vmifucalib and vmifustandard. The products 
of these were fed into vmifuscience which extracted the bias- 
subtracted, wavelength- and flux-calibrated, and relative fi- 
bre transmission-corrected science spectra. The final cube 
construction from the quadrant-specific science spectra out- 
put by the pipeline utilised an IFU table that lists the one to 
one correspondence between fibre positions on the IFU head 
and stacked spectra on the IFU (Bastian et al. 112006). A fi- 
nal data cube for each of the HRblue and HRorange grisms 
was then created by averaging over the final flux-calibrated 
science cubes created from each observation block. For a 
more detailed description of the data reduction and process- 
ing from science spectra to spatially reconstructed cube, see 
J09. A schematic representation of the data reduction pro- 
cesses performed within GASGANO is provided by Zanichclli 
et al. (20051. Sky subtraction was performed by locating a 



The cubes were corrected using an iRAF-based programme 
written by J. R. Walsh (based on an algorithm described in 



background region within each quadrant, summing the spec- 
tra over each of the spaxels in the reconstructed cube and 
subtracting the median sky spectrum of the region from its 
corresponding quadrant. Median combining was needed to 
ensure that any residual contamination from faint objects 
was removed. In the case of UM 420, the western edge of 
the IFU frame sampled the outermost faint regions of the 



unrelated foreground spiral galaxy UGC 1809 (e.g. Lopez 



Sanchez & Esteban (20081), and this background subtrac- 



tion procedure enabled the removal of any contaminating 
signal. Before averaging to produce a final data cube for 
each grism, each cube was corrected for differential atmo- 
spheric refraction (DAR). This correction accounts for the 
refraction of an object's spectrum along the parallactic angle 
of the observation as it is observed through the atmosphere. 



http: / / www.eso.org/ sci / data- processing/software / gasgano 



Walsh & Roy (1990)). The procedure calculates fractional 
spaxel shifts for each monochromatic slice of the cube rel- 
ative to a fiducial wavelength (i.e. a strong emission line), 
shifts each slice with respect to the orientation of the slit 
on the sky and the parallactic angle and recombines the 
DAR-corrected data cube. For both galaxies a high S/N ra- 
tio fiducial emission line (He I A5876) was available within 
the wavelength overlap of the HRblue and HRorange cubes. 
Although not essential, correcting multiple spectral cubes to 
a common emission line allows one to check that the align- 
ment between them is correct. 



2.3 Emission Line Profile Fitting 

Taking into account that there are 1600 spectra across the 
FoV of each data cube, we utilised an automated fitting pro- 
cess called PAN (Peak ANalysis; |Dimeo|[2005| . This is an 
IDL-based general curve-fitting package, adapted by |West^| 
moquette et al. (20071 for use with FITS data. The user 



can interactively specify the initial parameters of a spectral 
line fit (continuum level, line peak flux, centroid and width) 
and allow PAN to sequentially process each spectrum, fitting 
Gaussian profiles accordingly. The output consists of the fit 
parameters for the continuum and each spectral line's profile 
and the \ 2 value for the fit. It was found that all emission 
lines in the spectra of both UM 420 and UM 462 could be 
optimally fitted with a single narrow Gaussian. For high S/N 
ratio cases, such as Balmer lines and [O ill] A5007, attempts 
were made to fit an additional broad Gaussian profile under- 
neath the narrow component. These fits were assessed using 
the statistical F-test procedure, a process that determines 
the optimum number of Gaussians required to fit each line 
(see J09 for a more detailed discussion). We found that it 
was not statistically significant to fit anything more than 
a single Gaussian to the emission lines in the spectra of 
either object. Thus single Gaussian profiles were fitted to 
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each emission line profile, restricting the minimum FWHM 
to be the instrumental width. Suitable wavelength limits 
were defined for each emission line and continuum level fit. 
Further constraints were applied when fitting the [S n] dou- 
blet: the wavelength difference between the two lines was 
taken to be equal to the redshifted laboratory value when 
fitting the velocity component, and their FWHM were set 
equal to one another. The fitting errors reported by PAN 
underestimate the true uncertainties. We thus follow the er- 
ror estimation procedure outlined by J09, which involves the 
visual re-inspection of the line profile plus fit after checking 
which solution was selected by our tests, and taking into ac- 
count the S/N ratio of the spectrum. By comparing PAN fits 
to those performed by other fitting techniques (e.g. iraf's 
SPLOT task) on line profiles with an established configura- 
tion (i.e. after performing the F-test) we find that estimated 
uncertainties of ~5-10 per cent are associated with the sin- 
gle Gaussian fits. A listing of the measured flux and FWHM 
uncertainties appears in Table [3J where errors are quoted for 
individual component fits to emission lines detected on the 
integrated spectra from across all of UM 420 and UM 462. 



3 MAPPING LINE FLUXES 

3.1 Line fluxes and reddening corrections 

Full HRblue and HRorange spectra for UM 420 and UM 462 
are shown in Figure [TJa) and (b), respectively. Table [3] lists 
the measured FWHMs and observed and de-reddened fluxes 
of the detected emission lines within UM 420 and UM 462, 
respectively. The listed fluxes are from IFU spectra summed 
over each galaxy and are quoted relative to F(H/3) = 100.0. 
Foreground Milky Way reddening values of E(B — V) — 
0.04 and 0.02 were adopted, from the maps of Schlcgcl ct al. 
( 1998), corresponding to c(H/3) = 0.05 and 0.03 in the direc- 
tions of UM 420 and UM 462, respectively. The line fluxes 
were then corrected for extinction using the Galactic red- 



dening law of Howarth (19831 with Rv = 3.1 using c(H/3) 



values derived from each galaxy's Ha/H/3 and H7/H/3 emis- 
sion line ratios, weighted in a 3 : 1 ratio, respectively, after 
comparison with the theoretical Case B ratios from [Hum- 
mer fc Storey| ( |1987[ ) of jHa/ja/3 = 2.863 and j H7 /j H /3 = 
0.468 (at T e = 1CT K and N e = 100 cm" 3 ). The same 
method was used to create c(H/3) maps for each galaxy, us- 
ing Ha/H/3 and H7/H/3 ratio maps, which were employed 
to deredden other lines on a spaxel by spaxel basis prior to 
creating electron temperature, electron density and abun- 
dance maps. As stated in Table [3] overall c(H/3) values of 
0.25±0.08 and 0.19±0.04 were found to be applicable to the 
integrated emission from UM 420 and UM 462, respectively. 



3.2 Ha line properties 

Here we discuss the appearance of the targets in the light of 
the highest S/N ratio Hi emission line, Ha. Figure [2] shows 
flux, radial velocity and FWHM maps for Ha from UM 420 
and UM 462, respectively. 



3.2.1 UM 420 

The Ha flux map shown in the left-hand panel of Figure[5|a) 
displays a double peak; the first peak is centrally located at 
the RA and DEC listed in Table [T] with a surface brightness 
of 2.75x 10 -14 erg cm -2 s _1 arcsec -2 and the second peak is 
located ~2.2" east of that central position, with a peak sur- 
face brightness of 9.72xl0 -15 erg cm -2 s -1 arcsec -2 . These 
two Ha peaks provide the basis for defining two main star- 
forming regions, as displayed in Figure [3] An arm-like struc- 
ture can be seen north-west of the central peak, extending in 
a south-westerly direction; it could be the result of an ongo- 
ing merger or interaction between the two identified regions. 
A more luminous spiral galaxy, UGC 1809, lies 16.5" west- 
ward of the galaxy, as noted by |Takase ST Miyauchi -Isobej 
(19861), but at less than half the redshift of UM 420 it can- 



not be identified as an interacting companion. The VIMOS 
IFU Ha flux is a factor of 1.1 larger than the one measured 
by IT98 in a 3" x 200" slit and a factor of 1.8 times larger 
than that estimated by Lopez-Sanchez & Esteban ( 2008 \ 
from broad band Ha + [N 11] images. The Ha radial velocity 
map shown in the central panel of Figure [2] shows two sepa- 
rate velocity gradients. Firstly, there is a balanced velocity 
gradient at position angle (PA) 50° east of north (labelled 
as 'Position 1' on the central panel of Figure]^. This linear 
velocity gradient may be indicative of solid body rotation 
between the main emission peak (Region 1) and the 'arm- 
like' structure north-west of Region 1. Secondly, the main 
body of the galaxy, i.e. the double peaked emission struc- 
ture, shows a velocity gradient at a PA of 110° east of north 
(labelled as 'Position 2' on the central panel of Figure 

Figure [4] shows position- velocity (P-V) diagrams along 
both these axes. These diagrams were created by placing a 
spaxel-wide (i.e. with a width of 0.33") pseudo-slit across 
the radial velocity map, and plotting the radial velocity as a 
function of distance from the slit's centre. The smooth and 
symmetrically structured velocity gradient along Position 1 
suggests that the axis orthogonal to it could be adopted as 
a rotational axis for UM 420 (projected on the plane of the 
sky). From the top panel of Figure|4]we can then make a new 
estimate of the heliocentric systemic velocity of UM 420; a 
radial velocity of +17,545±30 km s -1 is needed to normalise 
the distribution along the proposed axis of rotation to zero 
(cf. the +17,514±12 km s" 1 estimation of IT98). Thus a ve- 
locity gradient ranging from ~-100 km s" 1 to ~+100 km s" 1 
is seen in the direction normal to the rotation axis, with the 
negative radial velocity of —100 km s -1 lying within the 
south-west arm of the galaxy. Also, a velocity gradient ex- 
ists between Region 1 and 2, ranging from ~ —5 km s -1 to 
~+100 km s . This velocity difference of ~100 km s -1 be- 
tween the two peaks in Ha flux may be evidence that Region 
2 is a satellite galaxy falling into or merging with a larger 
companion, Region 1. 

Another scenario could also be considered. Based on 
the Ha emission line map (Figure [3]), UM 420 could also be 
identified as a spiral galaxy with weak spiral arms (of which 
we detect only one arm - the south-west arm), a nuclear 
starburst (Region 1) and a faint, outer-lying starburst or 
superstar cluster situated at the edge of the galaxy (Region 
2). However, for a distance of 238 Mpc (Section [lj, UM 420 
has an outer radius of only ~3.5 kpc (at isophotes at 3% of 
the central peak flux), which would be very small for a spi- 
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Figure 1. VIMOS IFU summed spectra of UM 420 after smoothing with a 5 pixel boxcar function: (a) UM 420, top panel: high resolution 
blue spectrum, bottom panel: high resolution orange spectrum. The spectra are integrated over an area of 7.6"x6.3" (8.8x7.2 kpc 2 ) and 
correspond to exposure times of 4x402 s and 4x372 s for the HRblue and HRorangc grisms, respectively; (b) Same as (a), for UM 462. 
Spectra are integrated over an area of 10.2"xl0.2" (714.1x714.1 pc 2 ). 
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Table 3. Emission line measurements for summed spectra of UM 420 and UM 462. Observed and de-reddened fluxes are relative to 
F(H/3)=I(H/3)=100. Line fluxes were extinction-corrected using the c(H/3) values shown at the bottom of the table, calculated from the 
relative Ha, H/3 and H7 fluxes. FWHMs have been corrected for the instrumental FWHM resolution (Section \2.1\ 



Line ID 




FHWM (km sr 
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!) F(A) 
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!) F(A) 


1(A) 
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6548 [N II] 




238.3± 56.6 


10.9± 2.1 


9.1± 1.7 


127.3± 13.5 


3.1± 0.3 


2.7± 0.2 


6563 Ho 




178.8± 1.1 


334.4± 3.6 


279.1± 1.4 


103. 1± 0.6 


323.4± 2.5 


282. 2± 1.6 


6584 [N II] 




189.7± 14.4 


29.5± 1.8 


24.6± 1.5 


109.8± 4.2 


8.5± 0.3 


7.4± 0.2 


6678 He I 




135.7±105.2 


4.7±1.6 


3.8±1.3 


127.6± 6.3 


4.2± 0.4 


3.7± 0.3 


6716 [S II] 




196.5± 9.3 


33.9± 1.4 


27.9± 1.1 


101. 0± 5.2 


19.0± 0.8 


16. 4± 0.7 


6731 [S II] 




201.2± 12.6 


26.9± 1.4 


22. 1± 1.1 


105.9± 7.4 


14.4± 0.8 


12. 4± 0.7 


7136 [Ar III] 










104.6± 8.9 


9.2± 0.6 


7.8± 0.5 


c(H/3) dcx 






0.25±0.08 






0.19±0.04 




F(H/3)xl0" 16 erg s" 1 cm" 2 




359.1±6.1 






2660±14 





ral galaxy. Further to this, we find a spiral galaxy scenario 
improbable when considering the radial velocity map (cen- 
tral panel of Figure which clearly shows that Region 2 
is kinematically decoupled from the rest of the galaxy. 

A FWHM map for Ha is shown in the right-hand panel 
of Figure ^ a). The highest measured FWHM is ~270 km 
s~ , located southwards of the central emission peak, along 
the southern edge of the galaxy. The central flux peak is 
aligned with a peak of ~160 km s~ x whereas the second, 
easterly, flux peak is located within an area of decreased 
FWHM, ~130 km s — , surrounded by a slightly higher 
FWHM area of ~150 km s _ . The south-western arm har- 
bours the lowest line FWHM of ~80-120 km s" 1 . 



3.2.2 UM 462 

UM 462 exhibits a far more disrupted morphology than 
UM 420. The left-hand panel of Figure [2^b) shows a peak 
in Ha flux located at a = ll h 53 m 20.25 ± 0.1 s , S = 
-02° 32' 19.62 ± 0.1" (J2000), with a surface brightness of 
5.34xl0~ 14 erg cm -2 s" 1 arcsec -2 . Three additional flux 
peaks can be seen, located 2.5" northeast, 3.05" east and 
8.9" east of this central peak. These Ha peaks provide the 
basis for defining four main star-forming regions, as dis- 
played in Figure [5] Three of the regions have a similar peak 
surface brightness of 2.0-2.4X 10 -13 erg cm -2 s _1 arcsec -2 . 
Region 4 lies farthest away from the brightest SF region 



(Region 1), at almost 9"; a large low surface brightness en- 



velope (~5xl0 14 erg cm 2 s 



arcsec ) separates it from 



Regions 1-3. It could be the result of interaction with the 
nearby galaxy UM 461 (Taylor et al. 1995[ ). Previous ii-band 
and H I imaging by Taylor et al. ( 1995[ ) and VRI surface 
photometry by Tcllcs (1995) revealed only two star-forming 
regions within UM 462 which correspond to Regions 1 and 
2 in Figure [5] The Ha flux measured by IT98 amounts to 
35% of the corresponding VIMOS IFU flux, in accordance 
with the fraction of the galaxy's emitting area intercepted 
by IT98's slit. The radial velocity distribution of Ha, dis- 
played in the central panel of Figure [2jb) , shows no overall 
velocity structure through the galaxy. UM 462, being part of 
an interacting binary pair of galaxies, has a highly disturbed 
velocity distribution which shows no spatial correlation with 
flux (Figure[2jb): left-hand panel). The brightest region (Re- 
gion 1) has a radial velocity of ~0-10 km s~ , as does the 
majority of the emitting gas. The largest radial velocity lies 
directly south-east of this region, with a magnitude of +50 
km s _1 , in an area of low surface brightness emission. This 
velocity space extends into the envelope separating Region 
4 from Regions 1-3. Region 4 itself lies within a space of 
negative radial velocity, where again the peaks in velocity 
and flux are spatially uncorrelated. 

The Ha FWHM map of UM 462 (right-hand panel of 
Figure [5|b)) shows no spatial correlation with the flux dis- 
tribution and only a minimal correlation with radial velocity 
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Figure 2. Maps of UM 420 (a) and UM 462 (b) in Ha; left column: flux per 0.33 arcsec 2 spaxel; centre column: radial velocity (relative 
to heliocentric systemic velocity of +17,545 and +1057 km s — 1 for UM 420 and UM 462, respectively; right column: Ha FWHM corrected 
for the instrumental point spread function. See text for details. 



(Figure pb):central panel). A peak in the FWHM of Ha is 
•~ ~1 south of the emission peak of Region 4, at 130 km 
which aligns with the negative radial velocity peak. 



seen 
„— l 



distribution of [N n] is probably due to the effect of lower 
energy photons surviving farther away from the massive ion- 
ising clusters and allowed to produce N + beyond the 2+ 
zone. 



3.3 [O in] and [N n] 

Emission line maps of both galaxies in the light of [O ill] 
A5007 and [N n] A6584 lines are shown in Figs [6] and [7] For 
UM 420 the morphology of [O ill] is similar to that of Ha 
exhibiting two main starbursting regions along with the arc- 
like 'arm' protruding from the brighter region. In UM 462, 
which is the nearest of the two objects, the [O ill] map shows 
more structure that the corresponding Ha map with an ad- 
ditional area of compact emission to the south-west of Re- 
gion 1 of Fig [5] According to the morphological criteria of 
Cairos et al. (2001) these maps support the designation of 
UM 462 as being of 'iI,C type (irregular; cometary appear- 
ance) rather than of 4E' type; UM 420 can be classified as 
'il, M' (merging). For both BCGs, the [N 11] morphology is 
more diffuse and extended that that of [O ill]. The wider 



4 ELECTRON TEMPERATURE AND 
DENSITY DIAGNOSTICS 

The dereddened [O III] (A5007 + A4959)/A4363 intensity ra- 
tios were used to determine electron temperatures within 
each galaxy. The T c values were computed by inputting 
[O ill] intensity ratios and an adopted electron density into 
iraf'^TEMDEN task in the NEBULA package. Atomic tran- 
sition probabilities and collisional strengths for 2+ were 



taken from Wiese et al. ( 1996) and Lennon & Burke ( 19941, 



3 IRAF is distributed by the National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities for 
Research in Astronomy 
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Position 1 



Region 2 ^ Region 1 




-4 



-2 
arcsec 




Figure 3. UM 420: Ha emission per 0.33x0.33 arcsec 2 spaxel. 
The two main areas of Ha emission have been used to define two 
main star- forming regions. Also overlaid on the map are the loca- 
tions of position-velocity cuts used for the P-V diagrams shown 
in Figure [4] 



respectively. From the summed spectra over each galaxy 
(Table [3| and adopting the electron densities quoted be- 
low we derived mean electron temperatures of ~14,000 K 
and ~13,700 K for UM 420 and UM 462, respectively (Ta- 
ble]^. Using the same method, T c maps were derived using 
[O ill] (A5007 + A4959)/A4363 emission map ratios. Aver- 
age T c values for the regions corresponding to peaks in Ha 
emission for UM 420 and UM 462 are given in Tables |6] and 
[8] respectively, using the regionally integrated relative line 
intensities (Tables [5] The T c values from the summed 
spectra were adopted when computing jV c values from the 
[S 11] doublet ratio A6717/A6731: we derived average N e val- 
ues of ~170 cm -3 and ~90 cm -3 for UM 420 and UM 462, 
respectively (Table Q. Similarly, the T c maps were used to 
compute N B maps from A6717/A6731 ratio maps, which were 
then adopted for the computation of ionic abundance maps 
discussed in the following Section. Average N a values for the 
regions corresponding to peaks in Hcv emission for UM 420 
and UM 462 are given in Tables [6] and [8j respectively. 




-1 1 

Distance from FoV Centre / arcsec 




Radial Velocity km s 1 




-3 -2-10 1 

Distance from FoV Centre / arcsec 



5 CHEMICAL ABUNDANCES 

Abundance maps relative to H + were created for the N + , 
2+ , S + and S 2+ ions, using the AA6584, 5007, 6717+6731 
and 6312 lines, respectively, when the appropriate line was 
detected. Examples of flux maps used in the derivations of 
these abundance maps can be seen in Figures [6] and [7] show- 
ing the distribution of flux in the [O ill] A5007 and [N n] 
A6584 emission lines for UM 420 and UM 462, respectively. 
Since the [S in] A6312 line was not detected from UM 420, 
its S 2+ /H + abundance was estimated using an empirical re- 



Figure 4. Position velocity (P-V) diagrams for UM 420 in Ha: 
central panel: radial velocity map of UM 420 (see Figure |2ja) for 
details) where the black dotted line overlaid on the map represents 
cuts orthogonal to the two proposed axes of rotation, as used for 
the P-V diagrams in the top and bottom panels; top panel: P-V 
diagram for Ha along Position 1, defined as 50° east of north; 
bottom-panel: P-V diagram for Ha along Position 2, at 110° east 
of north. Zero-point velocities are relative to a recession velocity 
of +17,514 km s" 1 . Typical relative radial velocity errors are 
estimated as 30 km s _1 , with an estimated spatial uncertainty of 
0.5 spaxels (±0.17"). 
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Figure 6. UM 420: The distribution of flux per 0.33x0.33 arsec 2 spaxels of [O III] A5007 (left panel) and [N II] A6584 (right panel) 
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Figure 7. UM 462: The distribution of flux per 0.33x0.33 arsec 2 spaxels of [O III] A5007 (left panel) and [N II] A6584 (right panel) 



lationship between the S 2+ and S + ionic fractions (equation 
A38) from Kingsburgh fc Barlow] ( |1994[ ). Ionic abundances 
were calculated using iraf's IONIC task, using each galaxy's 
respective T e and -/V e maps described above, with each VI- 
MOS spaxel treated as a distinct 'nebular zone' with its own 
set of physical conditions. Ionic abundance ratios were also 
derived from the summed spectra fluxes, whose T 's and 
iV 's are listed in Table |4j using the program EQUIB06 (orig- 
inally written by I. D. Howarth and S. Adams). Table|4]also 
includes a Ne 2+ /H + abundance ratio for UM 420, as the 
[Ne ill] A3969 line was redshifted into the HRblue spectral 
range, and an Ar 2+ /H + ratio for UM 462. Abundance maps 
were however not created for either ion due to the low S/N 



ratio [Ne ill] A3969 and [Ar ill] A7135 emission line maps, 
and the N/H and Ar/H listed in Table[4]have been computed 
from integrated line ratios across the IFU. Differences be- 
tween abundances derived using the different atomic data 
tables utilised by IONIC and EQUIB06 were investigated by 
also using IONIC to compute abundances from the summed 
spectra. Typical differences were found to be ~2% or less. 
Since [O il] A3727 is not redshifted into the VIMOS spectral 
range for either galaxy we have to rely on published fluxes. 
We therefore adopted F(A3727)/F(H / 3) values of 2.02±0.04 
and 1.46±0.01 for UM 420 and UM 462, respectively (from 
IT98), and corrected them for reddening using the mean 
c(H/3) values derived here (see Table[3|. The spectra of IT98 
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Figure 5. UM 462: Ha emission per 0.33x0.33 arcsec 2 spaxel. 
The four main Ha emission areas have been labelled Regions 1—4. 



were obtained using the 2.1-m Kitt Peak Observatory Gold- 
Cam spectrograph's 3" x 200" slit aligned to collect the max- 
imum amount of flux from the targets. For UM 420, whose 
GoldCam and VIMOS Ha fluxes agree within 11 per cent, 
the adoption of [O n] fluxes from IT98 should introduce only 
minor uncertainties in our analysis; for UM 462 whose Ha 
flux measured by IT98 is 35% of the VIMOS flux, larger 
uncertainties might be expected. It is encouraging however 
that the relative line fluxes and intensities listed in Table [3] 
differ from those of IT98 by an average factor of only 0.96 
for UM 420 and 1.09 for UM 462. Hence for the derivation of 
the + abundance we adopted a constant J(A3727)//(H/3) 
ratio for each galaxy. As a result we cannot, for instance, 
comment on the spatial variance of the N/O abundance ra- 
tio across the galaxies. An alternative method for deriving an 
+ abundance map was investigated in which global + /N + 
ratios of 16.49 and 37.64 from IT98 for UM 420 and UM 462, 
respectively, were multiplied by the corresponding N + /H + 
IFU maps. Overall it was found that maps derived using 
the 'global factor' method yielded elemental oxygen abun- 
dances ~40% higher and nitrogen abundances ~55% lower 
than those derived from a constant /(A3727) //(H/3) ratio for 
each galaxy. However, by deriving an + abundance map in 
this way we would be prevented from creating a meaningful 
N/H map since ICF(N) = 0/0 + . We therefore opted to use 
the + abundance maps derived from the constant A3727 
fluxes when creating ICF(N), ICF(S) and O/H abundance 
maps. 

Ionic nitrogen, neon and sulphur abundances were con- 
verted into N/H, Ne/H and S/H abundances using ionisa- 



Table 4. Ionic and elemental abundances for UM 420 and UM 462 
derived, from summed IFU spectra over each galaxy. All ionisa- 
tion correction factors (ICFs) are taken from Kingsburgh & Bar- 
low, (1994). 



Property 



UM 420 



UM 462 



T e (O III)/ K 
N e (S II)/ cm" 3 
c(H3) 

He+/H+ (A4471) xlO 2 
He+/H+ (A5876) xlO 2 
He+/H+ (A6678) xlO 2 
He+/H+ mean XlO 2 
He 2 +/H+ (A4686)xl0 3 
He/H xlO 2 

0+/H+" xlO 5 
O++/H+ xlO 5 
O/H xlO 4 
12+log(0/H) 

N+/H+ xl0 e 
ICF(N) 
N/H xlO 6 
12+log(N/H) 
log(N/0) 

Ne++/H+ xlO 5 
ICF(Ne) 
Ne/H xlO 5 
12+log(Ne/H) 
log(Ne/0) 

S+/H+ xlO 7 
S++/H+ xlO 6 
ICF(S) 
S/H xlO 6 
12+log(S/H) 
log(S/0) 

Ar++/H+ xlO 7 
ICF(Ar) 
Ar/H xlO 7 
log(Ar/0) 



14000± 1500 
170± 80 
0.25±0.08 

7.65±1.03 
9.07±0.54 
10.6±3.6 
9.09±1.20 

9.09±1.20 

k u _i_3.35 
°- oi=t 1.56 

5 - 24± 113 

i-08±8:§| 
8.03±8i5 

1.94±H! 

i.95± 2 : J 2 

3.79±|;|| 
6.58± ; 44 

_l 4c_i_6.57 



13700± 350 

90± 75 
0.19±0.04 

6.68 ±0.62 
8.63 ±0.71 
10.4 ± 0.72 
8.58 ±0.72 
1.45 ± 0.08 
8.73 ±0.72 
4 qo + 0.43 
6.30±|]; 5 2 



1.06± 



().<)') 
I). OS 



n fii+ - 05 

u - oi=t 0.05 
Z.4Dzt 39 

l-50±8;fl 

6.i8±g:« 
-i-85±8:if 



3.33± 



2.24 
1.0!) 



2-05±J;i 
6.84±i 4 3 |8 
7-83±8;Sg 



-0.20± 



11.03 
0.01 



5-58±i'i2 
2-55 6 ±J;I| 

o.si±8:ig 

z - JO=c 1.84 
« 4Q_|_0.31 

_1 ficj+O- 43 
i - o,3=t 0.74 



'3 ijrv_|_0.33 

1.20±p| 
0.62±8;J| 
0.96±jj;| 
5.98±j,; 4 () 
-2.05±0; 22 

3.71±0 : }?s 

i-ot±8:8I 
3.98±8:lf 

-2.43±8;g? 



a Derived using [O II] A3727 flux from IT98 
b Derived using the relationship between S 2 + and S + (Equation 
A38) of Kinsgburgh & Barlow (1994) 



tion correction factors (ICFs) from Kingsburgh & Barlow 



(1994). Average abundances from integrated galaxy spec- 
tra are listed in Table [4] Elemental O /H abundance maps 
for UM 420 and UM 462 are shown in Figure ||J while av- 
erage ionic and elemental abundances derived from spectra 
summed over each individual star-forming region are listed 
in Tables ^ and [§] for UM 420 and UM 462, respectively^ 



Regional abundances derived directly by taking averages over 
the abundance maps were found to be inaccurate due to occasion- 
ally large T c uncertainties resulting from poor S/N ratio in some 
spaxcls of the [O III] A4363 maps. In contrast, summed regional 
spectra increased the S/N ratio for A4363, allowing more reliable 
average line ratios and hence temperatures and abundances for 
each region to be derived. 
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Table 5. UM 420 regional fluxes and intensities used for regional T e and N e diagnostics and regional ionic abundance calculations. 



Line ID 


Reg 


on 1 


Ref 


;ion 2 


F(A) 


1(A) 


F(A) 


1(A) 


4340 H7 


40.43± 0.49 


45.05± 0.54 


32.44± 1.15 


35.86± 1.27 


4363 [O III] 


4.51± 0.31 


5.00± 0.35 


5.34± 0.78 


5.88± 0.86 


4471 He 1 


4.14± 0.33 


4.49± 0.36 


4.92± 0.59 


5.30± 0.64 


4861 H/3 


100.00± 1.04 


100.00± 1.04 


100.00± 1.80 


100.00± 1.80 


4959 [O III] 


140.57± 1.47 


137.72± 1.44 


168.37± 2.69 


165.20± 2.64 


5007 [O III] 


441.01± 4.40 


427.64± 4.27 


514.34± 8.13 


499.91± 7.91 


5876 He 1 


13.04± 0.31 


10.86± 0.26 


11.29± 0.79 


9.53± 0.67 


6563 Ha 


366.84± 5.21 


278.98± 3.96 


296. 63± 3.88 


230.33± 3.01 


6584 [N n] 


29.66± 0.31 


22.50± 0.24 


16.80± 0.98 


13.02± 0.76 


6678 He 1 


2.24± 0.32 


1.68± 0.24 






6716 [S 11] 


30.15± 0.51 


22.52± 0.38 


24.11± 0.89 


18.40± 0.68 


6731 [S 11] 


23.31± 0.51 


17.37± 0.38 


20.81± 0.94 


15.86± 0.72 


F(H/3) a 


300.4± 3.330 


48.414± 0.874 



a in units of 10 16 erg s 1 cm 2 



Table 6. Ionic and elemental abundances for UM 420, derived 
from summed IFU spectra over each region. All ionisation correc- 
tion factors (ICFs) are taken from Kingsburgh & Barlow, (1994). 



Property 


Region 1 


Region 2 


T e (O III)/ K 


12300±330 


12900±680 


N e (S II)/ cm" 3 


120±55 


300±160 


c(H/3) 


0.37±0.04 


0.34±0.08 


0+/H+" xlO 4 


0.95± 0.12 


0.89± 0.25 


O++/H+ xlO 4 


0.79± 0.07 


0.92± 0.18 


O/H xlO 4 


1.74± 0.19 


1.81± 0.43 


12+log(0/H) 


8.24± 0.05 


8.26± 0.09 


N+/H+ xlO 6 


2.67± 0.18 


1.55± 0.22 


N/H xlO 6 


4.89± 1.60 


3.15± 2.39 


12+log(N/H) 


6.69± 0.12 


6.50± 0.25 


log(N/0) 


-1.55± 0.17 


-1.76± 0.33 


S+/H+ xlO 7 


13.38± 0.92 


11.26± 1.65 


S++/H+ 6 xlO 6 


7.47± 0.49 


6.39± 0.89 


S/H xlO 6 


7.58± 1.51 


5.87± 2.62 


12+log(S/H) 


6.88± 0.08 


6.77± 0.16 


log(S/0) 


-1.36± 0.12 


-1.49± 0.24 



a Derived using [O II] A3727 flux from IT98 
b Derived using the relationship between S 2 "*" and S + (Equation 
A38) of Kinsgburgh & Barlow (1994) 



5.1 UM 420 

The line fluxes for the summed spectra over Regions 1 and 
2 of UM 420 are listed in Table [5] and the average physical 
properties from which they are derived are listed in Table 
[6] A graphical representation of the variations in elemental 
and ionic abundances between Regions 1 and 2 is shown in 
Figure [9] 

The elemental oxygen abundance map shown in the 
left-hand panel of Figure [8] displays two peaks whose lo- 
cations correlate with the two Ha peaks of UM 420 (marked 
in Figure |3J. Minimal abundance variations are seen be- 
tween Regions 1 and 2, whose average metallicity is 12 + 
log(0/H) ~ 8.25 ± 0.07. Adopting a solar oxygen abun- 
dance of 8.71±0.10 relative to hydrogen ( |Scott et al.|2009[ ), 



this corresponds to an oxygen abundance of ~0.35 solar for 
Regions 1 and 2. This is ~0.2 dex higher than the value de- 
rived from the integrated spectrum of the whole galaxy. This 
could point to some limited oxygen enrichment in the areas 
corresponding to the peak Ha emission, that is the clusters 
associated with Regions 1 and 2. However, since all the oxy- 
gen abundances overlap within their la uncertainties, we 
cannot conclude that this variation in oxygen abundance is 
significant. The nitrogen abundances of the two regions of 
star-formation are consistent within their uncertainties (Ta- 
ble [6]). The N/O ratios for both regions are consistent with 
those of other metal-poor emission line galaxies of similar 
oxygen metallicity ( Izotov et al.|2006 1, and do not show the 
nitrogen excess reported by |Pustilnik et al.| ( |2004"| ). The S/O 
ratio for Region 1 is slightly higher than expected for BCGs, 
although Region 2 is closer to the average reported range of 
log(S/0) 1.4 - -1.7 (|Izotov et al."1[2006|) . 



5.2 UM 462 

Table[7]lists the fluxes from summed spectra over Regions 1- 
4 of UM 462. The average ionic and elemental abundances, 
along with the average T e and N e values, derived from the 
summed spectra over the respective regions are listed in Ta- 
ble [8] A graphical representation of the variation in elemen- 
tal and ionic abundances across Regions 1-4 is shown in 
Figure [9] 

The right-hand panel of Figure [8|a) shows the O/H 
abundance ratio map for UM 462. The oxygen abundance 
varies spatially across the different star-forming regions of 
UM 462 (as marked in Figure [5|. Four maxima in oxygen 
abundance are seen across the map, all aligning spatially 
with peaks in Ha emission (shown as overlaid contours), 
with the exception of Region 1, where the maximum in oxy- 
gen abundance appears to lie ~1.5" south of the Ha peak of 
Region 1. A decrease in the oxygen abundance can be seen 
~1.0" north-east of the peak in Region 1. Oxygen abun- 
dances for Regions 1-4 are listed in Table ||] A maximum 
variation of 40% is seen between Region 1 (displaying the 
highest metallicity) and Region 3 (displaying the lowest). 
The mean oxygen abundance of the four identified star- 
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Table 7. UM 462 regional fluxes and intensities used for regional T e and N e diagnostics and regional ionic abundance calculations. 



Region 1 Region 2 Region 3 Region 4 



Line ID 


F(A) 


1(A) 


F(A) 


1(A) 


F(A) 


1(A) 


F(A) 


1(A) 


4340 H7 


40.92± 0.72 


46.49± 0.49 


45.97± 0.41 


43.39± 0.54 


45.98± 0.41 


53.86± 0.48 


43.39± 0.54 


44.45± 0.56 


4363 [O in] 


8.61± 0.18 


9.26± 0.19 


7.27± 0.26 


7.55± 0.27 


9.17± 0.70 


10.67± 0.82 


6.62± 0.53 


6.77± 0.55 


4471 He 1 


3.01± 0.14 


3.19± 0.15 


3.51± 0.16 


3.61± 0.16 


3.66± 0.34 


4.12± 0.38 


3.63± 0.40 


3.70± 0.40 


4686 He n 






1.21± 0.12 


1.22± 0.12 










4861 H/3 


100.00± 0.65 


100.00± 0.65 


100.00± 0.65 


100.00± 0.65 


100.00± 0.65 


100.00± 0.65 


100.00± 0.76 


100.00± 0.76 


4959 [O in] 


218.90± 1.90 


215.76± 1.87 


169.17± 1.30 


167.95± 1.29 


182.95± 1.23 


177.54± 1.19 


174.63± 1.52 


173.83± 1.51 


5007 [O in] 


630.70± 5.95 


617.18± 5.82 


503.02± 2.64 


497.58± 2.61 


564.84± 2.32 


539. 98± 2.22 


507.96± 4.43 


504.49± 4.40 


5876 He 1 


14.47± 0.16 


12.72± 0.14 


9.77± 0.59 


9.16± 0.55 


15.02± 0.64 


11.48± 0.49 






6563 Ha 


332.75± 0.98 


274.41± 0.81 


325.84± 2.44 


295.81± 2.22 


483.79± 3.16 


324.20± 2.12 


293.19± 4.92 


275.84± 4.63 


6584 [N 11] 


4.25± 0.11 


3.50± 0.09 


6.88± 0.23 


6.24± 0.21 


9.88± 0.48 


6.60± 0.32 


5.04± 0.42 


4.74± 0.40 


6678 He 1 


2.34± 0.10 


1.91± 0.08 


3.21± 0.13 


2.90± 0.12 


5.19± 0.44 


3.41± 0.29 






6716 [S n] 


7.53± 0.10 


6.13± 0.08 


16.11± 0.20 


14.53± 0.18 


28.09± 0.72 


18.32± 0.47 


11.99± 0.38 


11.23± 0.35 


6731 [S II] 


5.56± 0.10 


4.52± 0.08 


11.91± 0.20 


10.74± 0.18 


19.90± 0.73 


12.95± 0.47 


8.79± 0.40 


8.23± 0.37 


F(H/3) a 


1065. 01± 6.32 


594.32±6.03 


96.24±2.33 


275.04± 15.62 



a in units of 10 16 erg s 1 cm 



Table 8. Ionic and elemental abundances for UM 462, derived from summed IFU spectra over each region. All ionisation correction 
factors (ICFs) are taken from Kingsburgh & Barlow, (1994). 



Property 


Region 1 


Region 2 


Region 3 


Region 4 


T e (O III)/ K 


13400±150 


13500±230 


15200±500 


12800±440 


N e (S II)/ cm" 3 


70±40 


70 ±40 




60±50 


c(H/3) 


0.26±0.03 


0.13±0.03 


0.54±0.03 


0.08±0.08 


0+/H+ a xlO 5 


4.85± 0.28 


4.38± 0.33 


3.91± 0.44 


5.09± 0.81 


0++/H+ xlO 5 


8.86± 0.30 


6.99± 0.34 


5.58± 0.49 


8.17± 0.88 


O/H xlO 4 


1.37± 0.06 


1.14± 0.07 


0.95± 0.09 


1.33± 1.49 


12+log(0/H) 


8.14-fc 0.02 


8.06± 0.02 


7.98± 0.04 


8.12± 0.05 


N+/H+X10 6 


0.34± 0.01 


0.60± 0.02 


0.50± 0.03 


0.51± 0.04 


N/H xlO 6 


0.97± 0.13 


1.56± 0.28 


1.22± 0.41 


1.33± 0.56 


12+log(N/H) 


5.99± 0.05 


6.19± 0.07 


6.08± 0.13 


6.12± 0.15 


log(N/0) 


-2.15± 0.07 


-1.86± 0.10 


-1.89± 0.17 


-2.00± 0.21 


S+/H+ xlO 7 


3.01± 0.09 


7.04± 0.28 


7.06± 0.45 


6.01± 0.51 


S++/H+ xlO 6 


0.96± 0.04 


1.27± 0.08 


1.37± 0.15 


1.97± 0.26 


S/H xlO 6 


0.66± 0.12 


1.15± 0.27 


1.32± 0.51 


1.49± 0.80 


12+log(S/H) 


5.82± 0.07 


6.06± 0.09 


6.12± 0.14 


6.17± 0.19 


log(S/0) 


-2.32± 0.09 


-2.00± 0.12 


-1.86± 0.19 


-1.95± 0.24 



a Derived using [O II] A3727 flux from IT98 



forming regions in UM 462 is 12 + log(0/H) = 8.08 ± 0.94, 
that is, ~0.25 solar. This is in good agreement with the value 
derived from the integrated spectrum of the whole galaxy. 
Regional N/H and S/H abundance ratios are given in Table 
|U Re gions 2, 3 and 4 have the same nitrogen abundance 
within the uncertainties. The maximum difference (0.2 dex) 
is between Regions 2 and 1; the latter shows the lowest nitro- 
gen abundance. In comparison with metal-poor emission line 
galaxies of similar metallicity ( Izotov et al.|[2006 l, UM 462 
appears to be overall slightly nitrogen poor. The sulphur 
abundance across the galaxy is variable with Region 1 hav- 
ing the lowest, a factor of ~ 2 below the other three regions. 
In comparison with other emission line galaxies, UM 462 dis- 
plays a lower than average log(S/0) of —2.03 [typical values 
range from ~ -1.4 to -1.7;|Izotov et al.| ([2006])]. 



5.2.1 Helium abundances 



Following the method outlined by |Tsamis et al.| ( |2003[ ) and 



Wesson et al. (20081, abundances for helium were derived 



using atomic data from Smits (19961, accounting for the ef- 



fects of collisional excitation using the formulae in Benjamin 
et al. (19991. Ionic and total helium abundances relative to 



hydrogen derived from helium recombination lines measured 
from summed spectra over the entire galaxies are given in 
Table [I] (derived using the global T c and N c values from 
the same Table). Given that the reddening varies across the 
galaxies (predominantly in the case of UM 462), and in or- 
der to minimise the uncertainties, we have also computed 
helium abundances per each individual star forming region; 
these are given in Tables [9] and [TTJ] for UM 420 and UM 462, 
respectively. These values are derived from summed spec- 
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Figure 8. O/H abundance ratio maps of UM 420 (left panel) and UM 462 (right panel) derived from their respective 2+ and + ionic 
abundance maps and T c and N c maps, as described in Section [5] Overlaid are the Ho flux contours from Figure [2fa). 



Table 9. Ionic and elemental helium abundances for UM 420 
derived from summed spectra over each region." 

Region 1 Region 2 



He+/H+ (A4471) xlO 2 


9.23±0.73 


10.77±1.23 


He+/H+ (A5876) xlO 2 


8.31±0.20 


7.15±0.52 


He+/H+ (A6678) xlO 2 


4.56±0.66 : 




He+/H+ mean XlO 2 


8.54±0.33 


8.05±0.70 


He 2 +/H+ (A4686)xl0 3 
He/H xlO 2 






8.54±0.33 


8.05±0.70 



a Entries followed by ':' were excluded from the final average. 



tra over each star-forming region of the galaxy dereddened 
with the corresponding c(H/3). Our adopted mean values 
for He+/H+ were derived from the A4471, A5876 and A6678 
lines, averaged with weights 1:3:1. The helium abundances 
were calculated for the average temperatures and densities 
given in Tables [6] and [8] Regional helium abundances were 
not calculated for Region 4 of UM 462 because only one He I 
line was detected (see Table [Tf. Helium abundances derived 
from summed spectra are found to agree with those derived 
in IT98, within the uncertainties. 



5.3 Stellar Properties 

5. 3. 1 Wolf-Ray et features ? 



Contrary to previous long slit studies ( Izotov & Thuan 1998 
|Schaerer et al.|1999||Guseva et al.|2000| ), which had limited 
spatial coverage, the VIMOS IFU spectra do not reveal any 
evidence for broad WR emission features at either ~4686 A 
or at ~5808, 5812 A. This is true for single spaxel spectra, 
summed spectra over Region 1 (the highest surface bright- 



ness SF region in both galaxies) or in the summed spectra 
over the whole of each galaxy, as shown in Figure flO] Hence 
this analysis does not support the classification of UM 420 
and UM 462 as Wolf-Rayet galaxies. When searching for 
WR features one must consider the width of the extraction 
aperture, which can sometimes be too large and thus dilute 



weak WR features by the continuum flux ( Lopez-Sanchez 
& Esteban 2008" I. We therefore also examined the summed 



spectra over each individual star-forming region, with ex- 
traction apertures equal to those displayed in Figures [3] and 
[5] aiming to reduce the dilution of any WR features. How- 
ever, those remained undetected. This supports the state- 



ment by Schaerer et al. ( 1999) regarding the dependence of a 



object's classification as a 'WR galaxy' on the quality of the 
spectrum, and location and size of the aperture. Given that 
our spectroscopy was obtained with an 8.2 m telescope, with 
exposure times of ~1500 s for each spectral region, we con- 
clude that the presence of WR stars in UM 420 and UM 462 
is not yet proven. 



5.3.2 Starburst ages and star formation rates 

Luminosities of hydrogen recombination lines, particularly 
H/3, can provide estimates of the ionising flux present, 
assuming a radiation-bounded nebula ( |Schaerer fc Vacca| 
|1998[ ). Thus, the equivalent width (EW) of H/3 is commonly 
used as an age indicator of the ionising stellar population 
at a given metallicity. Maps of MU(H/3) for UM 420 and 

: (b), respectively. 



UM 462 can be seen in Figure 11a 
We can use these maps, in conjunction with the metallicity 
maps described in Section [5] to estimate the age of the lat- 
est star forming episodes throughout UM 420 and UM 462 
by comparing the regional observed average EW(R(5) values 
with those predicted by the spectral synthesis code STAR- 
BURST99 ( |Leitherer et al.||1999[ ). For the models we chose 
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Table 10. Ionic and elemental helium abundances for UM 462 derived from summed spectra over each Region (with the exclusion of 
Region 4, where only one He I line was detected)." 





Region 1 


Region 2 


Region 3 


He+/H+ (A4471) xlO 2 


6.66±0.32 


7.54±0.34 


8.75±0.81 


He+/H+ (A5876) xlO 2 


9.97±0.10 


7.20±0.43 


9.33± 0.40 


Hc+/H+ (A6678) xlO 2 


5.34±0.23 : 


8.11± 0.33 


9.83±0.84 


He+/H+ mean XlO 2 


9.14±0.16 


7.45±0.39 


9.31±0.57 


He 2 +/H+ (A4686)xl0 3 
He/H xlO 2 




1.04±0.11 




9.14±0.16 


7.55±0.40 


9.31±0.57 



Entries followed by ':' were excluded from the final average. 



UM420 All Regions 



UM420 All Regions 




UM462 All Regions 



UM462 All Regions 




4600 4700 

Rest Wavelength / A 



5800 

Rest Wovelength / A 



Figure 10. Sections of VIMOS IFU spectra, in flux units of X 10 — 16 erg s — 1 cm -2 , where the WR 'blue bump' feature (left column) and 
C IV 'red' feature (right column) would be located within summed spectra over the core region (labelled Region 1 in Figures [3] and [5] for 
UM 420 and UM 462, respectively) and over the entire galaxy. There are no detections of WR emission features in any of the spectra. 



metallicities of 0.4 and 0.2 Zq, representative of the average 
metallicities of UM 420 and UM 462 (derived from regional 
summed spectra), together with assumptions of an instan- 
taneous burst with a Salpeter initial mass function (I MF), 
a total mass of 1 X 10 6 Mq (the default mass chosen by Lei- 



therer et al. (19991 to produce properties that are typical 
for actual star- forming galaxies), and a 100M e upper stellar 



mass limit (which approximates the classical Salpeter ( 1955 1 



IMF). Models were run for Geneva tracks with "high" mass 
loss rates (tracks recommended by the Geneva group) and 
Padova tracks with thermally pulsing AGB stars included. 
For each of these evolutionary tracks, two types of model at- 
mosphere were used; firstly the Pauldrach-Hillier (PH, the 
recommended atmosphere) and secondly Lejeune-Schmutz 
(LS). The latter was chosen because it incorporates stars 
with strong winds, which would be representative of the WR 
population within each galaxy (if present). The stellar ages 
predicted by each model and the observed average £W(H/3) 
within each peak emission region in UM 420 and UM 462 
were derived from Figure [12] and are listed in Table [IT] The 
difference between the ages predicted by the Geneva and 



Padova stellar evolutionary tracks is relatively small, with 
the Geneva tracks predicting lower ages by up to 20%. The 
difference in ages predicted by the PH and LS atmospheres 
are smaller still, with LS predicting lower ages by ~6%. We 
cannot comment on which atmosphere model is more ap- 
propriate given that the existence of WR stars within each 
galaxy is questioned by the current study, hence we adopt 
average stellar ages from the four model combinations. The 
current (instantaneous) star formation rates (SFRs) based 
on the Ha luminosities were calculated following Ke nnicutt| 
(19981 and are given in Table 11 SFRs corrected for the 



sub-solar metallicities of these targets are also given, derived 



following the methods outlined in Lee et al. (20021. 

UM 4%0: Regions 1 and 2 contain ionising stellar 
populations with weighted average ages of 3.94±0.31 and 



4.99±0.21 Myr, respectively (Table 111. Interestingly, the 



ages returned by all 4 models consistently show Region 2 
to be older than Region 1, which could imply that these 
were in fact once isolated bodies, and that star formation 
within Region 1 was triggered or is still a product of an on- 
going merger with Region 2. This also agrees well with the 
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Table 11. Age of the latest star formation episode and current star formation rates. 
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Starburst ages (Myr) 



UM 420 



UM 462 



Model 

Padova-AGB PH 
Padova-AGB LS 
Geneva-High PH 
Geneva-High LS 



Reg. 1 
4.51±0.21 
4.24±0.34 
3.76±0.36 
3.50±0.34 



Reg. 2 
5.26±0.19 
5.16±0.19 
4.91±0.27 
4.71±0.22 



Reg. 1 
4.85±0.32 
4.66±0.36 
3.91±0.56 
3.76±0.48 



Reg. 2 
5.45±0.13 
5.31±Q.10 
4.86±0.13 
4.61±0.10 



Reg. 3 
5.85±0.13 
5.66±0.07 
5.21±0.10 
5.01±0.10 



Reg. 4 
5.90±0.27 
5.76±0.24 
5.51±0.39 
5.21±0.33 



Star formation rates (Mq yr 1 ) 



UM 420 



UM 462 



SFR(H«) a 
SFR(Ho) 6 



10.5 
7.0 



1.31 
0.87 



0.104 
0.069 



0.047 
0.027 



0.021 
0.012 



0.018 
0.013 



a Derived using the relationship between SFR and L(Ho) from 
b Corrected for sub-solar metallicities following Lee et al. ( 2002 I 



Kennicutt| | |1998| l. 



in-falling merger scenario derived from radial-velocity maps 



(discussed in Section 3.2.1 1. However, the significance of the 



derived difference in ages between the two regions is only 
at the 2a level. Also in support of this are their respective 
SFRs which differ by almost an order of magnitude. Global 
SFRs derived by |L6pez-Sanchez fc Este ban ( 2008J) are 3.7 
± 0.2 M© yr -1 (based on an Ha flux estimated from broad- 
band photometry) and 1.9 ± 0.9 Mq yr -1 (based on the 
1.4 GHz flux), and are lower than our measurements which 
are based on the monochromatic Ha flux. 

UM 462: Table [7T1 show average stellar population ages 
of 4.19±0.42, 5.07±0.11, 5.45±0.10 and 5.56±0.31 Myr for 
Regions 1-4, respectively. Since Region 1 shows the highest 
ionising flux, we would expect it to contain the youngest 
stellar population. It is notable that the age of the stellar 
population increases through Regions 2-4 as the distance of 
separation from Region 1 increases and at the same time the 
SFR decreases. Our SFRs are comparable to previously pub- 
lished values of 0.13 (based on FIR fluxe s), 0.00081 (blue' ) 
and 0.5 Mq yr - (narrow-band Ha) from Sage et aL (1992 i 
and . 36 (60 /im) and 0.13 M Q yr" 1 (1.4 GHz) from |Hopkins 
et al. (20021. Evidence suggests that Regions 1-3 are linked 



in their SF properties, with decreasing metallicity and SFR 
and increasing stellar population ages as the distance from 
Region 1 increases. Region 4 breaks this pattern in metal- 
licity harbouring a starburst which is coeval with that of 
Region 3. 



6 EMISSION LINE GALAXY 
CLASSIFICATION 

It is useful to attempt to diagnose the emission line excita- 
tion mechanisms in BCGs using the classic diagnostic dia- 
grams of |Baldwin et al.|fl981| > (the BPT diagrams). These 
are used to classify galaxies according to the dominant exci- 
tation mechanism of their emission lines, i.e. either photoion- 
isation by massive stars within H II regions or photoionisa- 
tion by non-thermal continua from AGN. The diagrams con- 
sist of excitation-dependent and extinction-independent line 
ratios: log([0 ill] A5007/H/3) versus either log([S II] A6716 
+ A6731)/Ha) or log([N n] A6584/Ha). Star-forming galax- 



ies fall into the lower left region of the diagrams, AGN host 
galaxies fall into the upper right region and Low-Ionisation 
Emission Line Regions (LINERs) fall on the lower right. 
Here we adopt the 'maximum starburst line' derived by 



Kewley et al. (20011 from starburst grids defined by two 
parameters; metallicities Z = ~ 0.05-3.0 Zq and ioniza- 
tion parameters q = 5 x 10 6 - 3 x 10 s cm s _1 , where q 
is the maximum velocity of an ionisation front that can be 
driven by the local radiation field and is related to the non- 
dimensional ionisation parameter U = q/c, ( |Dopit; 



et al. 



20011. This line defines the maximum flux ratio an object 



can have to be successfully fitted by photoionization mod- 
els alone. Ratios lying above this boundary are inferred to 
require additional sources of excitation, such as shocks or 
AGNs. It is thus instructive to see what area of the BPT 
diagram UM 420 and UM 462 occupy as a whole, as well 
as their resolved star-forming regions which were defined in 
Figures [3] and [5] 



UM 420: Figure 13 » and (b) shows the BPT diagram 



locations of the two star-forming regions within UM 420. Al- 
though a large spread can be seen for both regions, Region 
2 occupies the higher end of the [O ill] A5007/H/3 emis- 
sion line ratio in both diagrams, and does not exceed the 
[O ill] A5007/H/3 values of Region 1. Some spaxels do cross 
the 'maximum starburst line', but are insufficient in fre- 
quency to be considered as evidence for substantial non- 
thermal line excitation. The [S ll]/Ha ratios are generally 
lower than those predicted by standard predictions of shock 
models from the literature. There is a larger spread in [N n] 
A6584/Ha as compared to [S ll]/Ha which is probably due 
to the fact that the N + /H + abundance ratio differential 
between Regions 1 and 2 is larger than the corresponding 
S + /H + differential by 45 per cent. In relation to the grid of 



(Z, q) models of |Kewley et al.| ( |200l[ ), UM 420 lies within a 
range of Z = 0.4-0.5 Z Q and q > 3 x 10 8 . This is consistent 
with the average metallicity of 0.35 ± 0.06 Zq for Regions 1 
and 2 (Table IHJ. We conclude that photoionisation by stel- 
lar populations is the dominant line excitation mechanism 
within UM 420. 

UM 462: The BPT diagnostic diagrams shown in Fig- 
ure |14[ a) and (b) show the emission line ratios for each 
spaxel within the four star-forming regions of UM 462 (as 



16 B. L. James et al. 



UM 420 




N+ N 0+ 02+ O S+ S2+ S N/0 S/0 

UM 462 



o 

-t— ■ 

CO 

I 

CD 
C_> 

c 
cC 
"O 

c 
13 
n 
co 

I 



10" 
10 3 
10 A 
10 5 
10 e 
10 7 



■ 


Reg 1 


• 


Reg 2 


▲ 


Reg 3 


V 


Reg 4 



* 



N+ N 0+ 02+ O S+ S2+ S N/O S/O 



Figure 9. Graphical representation of the variations in ionic and 
elemental abundances across UM 420 and UM 462, as given in 
Tables [6] and [8] respectively. 
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defined in Figure [5). As with UM 420 the [S Il]/Ha diag- 
nostic ratios tend to straddle the theoretical upper limit 
for photoionisation. Although their locations appear highly 
clustered, a slight distinction can be seen for Region 4 which 
lies slightly to the left of the diagram. The large spread in 

b), where the 



14 



values can be seen more clearly in Figure 
spread along the X-axis, i.e. [N n] A6584/Ha, extends on ei- 
ther side of the other three regions. This range in [N n] and 
[S n] excitation conditions is perhaps to be expected given 
the 9" distance separating Region 4 from the central regions. 
Several spaxels' emission line ratios do spread into the locus 
of non-thermal excitation but their low occurrence renders 
them insignificant. In relation to the models of |Kewley et al.| 
(20011), UM 462 lies within a range of Z =0.2-0.5 Z and 



q > 1.5 x 10 . This is consistent with the average metallic- 
ity of ~O.24Z for Regions 1-4 (Table Isb. We conclude that 
photoionisation from stellar sources is the dominant exci- 
tation mechanism within UM 462 as well. Comparing the 
two galaxies, the data points of UM 462 are more highly 
clustered than those of UM 420, despite the former system's 
disrupted morphology and are further to the left (smaller 
values) on the X-axis of both diagrams. The explanation 
probably lies in the higher S/H and N/H abundances of 



Figure 11. Maps of the equivalent width of H/3 across: (a) 
UM 420; (b) UM 462. Overlaid red solid lines are contours of 
EW(Hf3); black-dashed lines are Ha flux contours from Figure [5] 



UM 420 compared to those of UM 462 (by corresponding 
factors of ~ 2.5). Both galaxies display a similar level of 
[O ill] A5007/H/3 excitation. 



7 CONCLUSIONS 

We have analyzed VIMOS IFU integral field spectroscopy 
of the unrelated BCGs UM 420 and UM 462 and stud- 
ied their morphology by creating monochromatic emission 
line maps. Both systems show signs of interaction and/or 
perturbation with the former galaxy currently undergoing a 
merger (type 'iI,M') and the latter displaying a highly dis- 
rupted irregular or cometary appearance (type 'iI,C; prob- 
ably due to interaction with a BCG which was not part of 
this study, UM 461). The spatially resolved emission line 
maps in the light of [O in] and Ha have revealed two main 
areas of massive star formation in UM 420 along with an 
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Figure 12. EW(Hf3) as a function of age, as predicted by the 
STARBURST99 code for metallicities of O.4Z (UM 420) and 
O.2Z0 (UM 462) using a combination of Geneva or Padova stel- 
lar evolutionary tracks and Lejeune-Schmutz (LS) or Pauldrach- 
Hillier (PH) model atmospheres. The observed average i?lU(H/3)'s 
for each star forming region in UM 420 (top panel) and UM 462 
(bottom panel) are overlaid (solid line). 



'arm-like' structure, and at least four such areas in UM 462. 
Current star formation rates were computed from the Ha 
line luminosities for each main starbursting region and the 
ages of the last major star formation episode were estimated 
by fitting the observed Balmer line equivalent widths with 
STARBURST99 models. The two merging components of 
UM 420 have SFRs that differ by a factor of ~ 8 and star- 
burst episodes separated by 1 Myr. The latest major star 
formation event took place ~ 4 Myr ago in the largest merg- 
ing component, and has been producing stars at a rate of 
10 Mq yr _1 . In UM 462 the last star forming episode was ~ 
4-5 Myr ago and has been producing stars at a rate which 
varies across the galaxy between ~ 0.01 - 0.10 Mq yr~\ 
indicative of propagating or triggered star formation. 

For both targets the abundances of He, N, O, and S 
were measured and O/H abundance ratio maps were created 
based on the direct method of estimating electron tempera- 
tures from the [O ill] A4363/A5007 line ratio. The measured 
oxygen abundances are 12 + log(0/H) = 8.03 ± 0.20 for 
UM 420 and 8.03 ± 0.10 for UM 462 (20 per cent solar). We 
find no evidence for significant nitrogen or oxygen variations 
across the galaxies (at the 0.2 dex level), which would point 
to self-enrichment from nucleosynthetic products associated 
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Figure 13. Emission line diagnostic diagrams for UM 420. Each 
data point represents a spaxel in the ratioed dereddened flux 
maps corresponding to [O in] A5007/H/3 versus (a) [S II] A6716 
+ A6731/Hq and (b) [N II] A6584/Ha. Spaxels corresponding to 
the two star- forming regions defined in Figure [3] are separated by 
colour. The empirical diagnostic 'maximum starburst' line from 
|Kewley et al.| ( |2001| l is also shown, as are the positions of emis- 
sion line ratios for star-forming galaxies(SFG) and active galaxies 
(AGN). 



with the recent massive star formation activity. Regarding 
the abundance of nitrogen, and the N/O ratio, this result is 
in qualitative agreement with the finding that these BCGs 
cannot be classified as Wolf-Rayet galaxies as the charac- 
teristic broad-line stellar features have not been detected by 
VIMOS; the existence of WR stars in them is an open issue. 
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UM 462 Emission Line Diagnostics 
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Figure 14. Same as for Figure |13| but for UM 462. Regions 
1-4 correspond to the separate star-forming regions denned in 
Figure [5] 
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